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Abstract 

We propose a scenario of the right-handed neutrino dark matter in the context of the minimal 
gauged U(1)b-l model by introducing an additional parity which ensures the stability of dark 
matter particle. The annihilation of this right-handed neutrino takes place dominantly through the 
s-channel Higgs boson exchange, so that this model can be called Higgs portal dark matter model. 
We show that the thermal relic abundance of the right-handed neutrino dark matter with help of 
Higgs resonance can match the observed dark matter abundance. In addition we estimate the cross 
section with nucleon and show that the next generation direct dark matter search experiments can 
explore this model. 
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I. INTRODUCTION 



The nonvanishing neutrino masses have been confirmed by various neutrino oscillation 
phenomena and indicate the evidence of new physics beyond the Standard Model. The most 
attractive idea to naturally explain the tiny neutrino masses is the seesaw mechanism in 
which the right-handed (RH) neutrinos singlet under the SM gauge group are introduced. 
The minimal gauged U(l) b-l model based on the gauge group SU(3)c x SU (2) L x U(l)y x 
U(1)b-l is an elegant and simple extension of the SM, in which the RH neutrinos of 
three generations are necessarily introduced because of the gauge and gravitational anomaly 
cancellations. In addition, the mass of RH neutrinos arises associated with the U(1)b-l 
gauge symmetry breaking. 

Although the scale of the B — L gauge symmetry breaking is basically arbitrary as long as 
phenomenological constraints are satisfied, one interesting option is to take it to be the TeV 
scale {3)]. It has been recently pointed out ^| that when the classical conformal invariance 
is imposed on the minimal U(1)b-l model, the symmetry breaking scale appears to be the 
TeV scale naturally. If this is the case, all new particles, the Z' gauge boson, the B — L 
Higgs boson H and the RH neutrinos appear at the TeV scale unless the U(1)b-l &^E e 
coupling is extremely small, and they can be discovered at Large Hadron Collider [5|-l8]. 
Then we may be able to understand the relation between the gauge symmetry breaking and 
the origin of neutrino masses. 

Although such a TeV scale model is interesting and appealing, one might think that the 
absence of dark matter (DM) candidate is a shortcoming of this model. A sterile RH neutrino 
with mass of the order of MeV is one possibility (9]. In this paper, we propose a very simple 
idea to introduce the DM candidate in the minimal gauged U(1)b-l model. We introduce 
the Z 2 parity into the model and impose one of three RH neutrinos to be odd, while the 
others even. In this way, the Z 2 -odd RH neutrino becomes stable and the DM candidate. 
Note that two RH neutrinos are enough to reconcile with the observed neutrino oscillation 
data, with a prediction of one massless light neutrino. Therefore, without introducing any 
additional new dynamical degrees of freedom, the DM particle arises in the minimal gauged 
U(1)b-l model. 

The paper is organized as follows. In the next section, we briefly describe our model. In 
section III, we estimate the thermal relic density of the RH neutrino and identify the model 
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parameter to be consistent with the current observations. We also calculate the scattering 
cross section between the DM particle and nucleon and discuss the implication for the direct 
DM search experiments. We summarize our results in the section IV. Our notations and the 
formulas used in our analysis are listed in Appendix. 



II. THE MINIMAL GAUGED U(1) B -l MODEL WITH Z 2 PARITY 

The model is based on the gauge group SU(3) C x SU(2) L x U(1) Y x U(1) B -l- Additional 
fields besides the standard model fields are a gauge field Z'^ of the U(1) B -l, a SM singlet 
B — L Higgs boson \1> with two U(1) B -l charge, and three RH neutrinos N { which are 
necessary for the gauge and gravitational anomaly cancellations. In describing the RH 
neutrinos, we use the four component representation of RH neutrino constructed from the 
Weyl spinor u Ri1 

*-(.*)• 

For the two RH neutrinos, Ni and N 2 , we assign Z 2 parity even, while odd for N 3 , so that 
the RH neutrino iV 3 is stable and, hence, the DM candidate. 

Due to the additional gauge symmetry U(1)b-l, the covariant derivative for each fields 
is given by 

D, = - iq B _ L g B _ L Z^ (2) 

where is the covariant derivative in the SM, and q B -L is the charge of each fields 

under the U(1) B -l with its gauge coupling g B -L- 

Yukawa interactions relevant for the neutrino masses are given by 

3 2 ^3 
C mt = J2J2 V<* l °& N i ~ 2 Yl X RA^ p R N i + h -c-, (3) 

a=l i=l i=l 

where $ = — ir 2 $* for $ being the SM Higgs doublet, and without loss of generality we have 
worked out in the basis where the second term in the right-hand-side is in flavor diagonal 
for RH neutrinos. Because of the Z 2 parity, the DM candidate A 3 has no Yukawa couplings 
with the left-handed lepton doublets. 

The general Higgs potential for the SU (2) l doublet $ and a singlet B — L Higgs ^ is 
generally given by 

V(Q,V) = m\\<5>\ 2 + m\\^\ 2 + + A 2 |^| 4 + A 3 |$| 2 |^| 2 . (4) 



The Higgs fields and ip are obtained by expanding $ and 9 as 

$ = ( ° ] , (5) 

9 = -L(v' + (6) 

around the true vacuum with the vacuum expectation values v and v'. These are related 
with the mass eigenstates h and H through 



(7) 



cos 9 — sin 9 
sin 9 cos 9 I \ ip 

with 9 being the mixing angle. Their masses are given by 

M 2 h = 2X 1 v 2 cos 2 9 + 2\ 2 v' 2 sm 2 9~2X 3 vv'sm9cos9, (8) 
M% = 2X 1 v 2 sm 2 9 + 2X 2 v' 2 cos 2 9 + 2X 3 vv'sm9cos9. (9) 

The mass of the new neutral gauge boson Z' arises by the U(1)b-l gauge symmetry 
breaking, 

M% = 4g 2 B _ L v' 2 . (10) 
Associated with the U(l) b-l gauge symmetry breaking, the RH neutrinos Ni acquire masses 

M Ni = -\ R ^. (11) 

p experiment, the current lower bound on the Z' boson mass has been found to 

2v' > 6 - 7 TeV. (12) 



From LE 
be 
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M z , 



9B-L 

Two ^2-even RH neutrinos N\ and N 2 are responsible for light neutrino masses via the 
seesaw mechanism, 

V 2 

i=l,2 Ni 

Note that the rank of this mass matrix is two, so that the lightest neutrino is massless. 



III. RIGHT-HANDED NEUTRINO DARK MATTER 



Due to the Z 2 parity, one of RH neutrino A^ 3 (we denote it as N hereafter) in our model 
can be the DM candidate. We first estimate its relic abundance and identify the model 

4 



parameters to be consistent with the current observations. Next we calculate the scattering 
cross section between the DM particle and a proton and discuss the implication for the direct 
DM search experiments. 



A. Thermal relic density 

The DM RH neutrino interacts with the SM particles through couplings with B — L 
gauge and B — L Higgs bosons. Note that neutrino Dirac Yukawa interactions are absent 
because of the Z 2 parity. The most of annihilation of the RH neutrinos occurs via Z', H and 
h exchange processes in the s-channel. In practice, the dominant contributions come from 
the Higgs (h and H) exchange diagrams, because the Z' exchange processes are suppressed 
by the inverse square of the B — L Higgs VEV v ' > 3 TeV. Thus, we obtain Higgs portal DM 
of RH neutrino effectively. The relevant annihilation modes are the annihilation into //, 
W + W~, ZZ, and h(H)h(H). Since RH neutrino DM couples to only B — L Higgs * while 
a SM particle does to SM Higgs $, the DM annihilation occurs only through the mixing 
between these two Higgs bosons. Although it is not so severe, the precision electroweak 
measurements [3] as well as the unitarity bound [3] give constraints on the mixing angle 
and mass spectrum of the Higgs bosons. 

The thermal relic abundance of DM 

n N h 2 = 1.1 x 10 9 m f/ T " . GeV 1 , (14) 

with the Planck mass Mp, the thermal averaged product of the annihilation cross section 
and the relative velocity (crv), the total number of relativistic degrees of freedom in the 
thermal bath g+, and the decoupling temperature T^, is evaluated by solving the Boltzmann 
equation for the number density of RH neutrino n^; 

( ^ + 3Hn N = -(av)(n 2 N -n 2 EQ ), (15) 

and the Friedmann equation 
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a \ 87T 



H ^{l)-3k P ' (16> 
with ?2eq and a(t) being the equilibrium number density and the scale factor, under the 
radiation dominated Universe with the energy density p = p ra a [yj . 



Fig. 1 shows the relic density Q^h 2 as a function of the DM mass for a set of 
parameters: (V, M h , M H , M z >, sin 6) = (4000 GeV, 120 GeV, 200 GeV, 1000 GeV,0.7), for 



example. Will 



dnson Microwave Anisotropy Probe measured the value of DM abundance as 



Vt DM h 



0.1 



15] . The figure shows that a desired DM relic abundance can be obtained for 



only near Higgs resonances, ~ M^/2 or M#/2. 

Fig. 2 shows the relic density Vl^h 2 as a function of the DM mass mjy for a smaller Higgs 
mixing sin# = 0.3 (others are the same as in Fig. 1). Compared with Fig. [TJ for < My/ 
where the DM particles dominantly annihilate into //, the relic density further increases 
because of the small mixing angle. When the DM is heavier, the annihilation mode into 
Higgs boson pairs is opened and the relic density slightly deceases, but the reduction is not 
enough to reach Q^h 2 ~ 0.1. 
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FIG. 1: The thermal relic density of RH neutrino DM as a function of its mass for a parameter 
set: (v',M h ,M H ,M Z ',sm9) = (3000 GeV, 120 GeV, 200 GeV, 1000 GeV, 0.7). 



Our model is quite analogous to the so-called gauge singlet scalar dark matter [16H18]. 



Some recent studies can be found in Refs. 
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20] . In the gauge singlet scalar DM model, the 



thermal abundance is mainly controlled by the interactions between the SM Higgs boson and 
the DM particle. In our model, B — L Higgs VEV v' can play the same role for tjin < My/, 
namely a larger v' corresponds to weaker coupling between DM and Higgs for a fixed DM 
mass. On the other hand, for > My/ the difference appears. Even if the annihilation 




FIG. 2: The same as Fig. Q] but for sin6» = 0.3. 

mode into VF-boson pair becomes kinematically available, it is not possible to obtain the 
desired DM abundance without the Higgs resonant annihilation because the bound on v' 
given by Eq. (|T2|) is stringent. 

B. Direct detection of dark matter 



2l| and studied [22|, |23|. The main process 



Our RH neutrino DM can elastically scatter off with nucleon, unlike another RH neutrino 
DM model has been proposed by Krauss et. al. 

is Higgs exchange and the resultant cross section for a proton is given by 

4 / m p m N \ 2 2 
si ~~ _ I ™ i ™ J J pi 



a 



7r \m p + m N 



with the hadronic matrix element 



f P _ f {p) a q 2 (p) \- a q 

m„ ' y m„ 27 ^ m 



q=u,d,s 



c,b,t 



and the effective vertex (see Appendix for notations) 



'<9$ 1 » 9$ 1 



where m g is a mass of a quark with a Yukawa coupling y 9 , and and are constants. 
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(17) 



(18) 



(19) 



From Eq. (119)). one can see that oc (sin2#/i/) 2 for a given DM mass m^- Fig. 3 shows 
the spin-independent cross section of RH neutrino with a proton. The resultant cross section 



is found to be far below the current limits reported by XENON10 [24] and CDMSII |25] : 
0"si ^ 4 x 10 -8 — 2 x 10~ 7 pb, for a DM mass of 100 GeV-1 TeV. Future experiments such 
as XENON1T 26] can reach the cross section predicted in our model. 
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FIG. 3: The spin independent scattering cross section with a proton. All parameters are same as 
those used in the previous section. The upper and lower lines correspond to sin# = 0.7 and 0.3, 
respectively. 



IV. SUMMARY 

We have proposed a scenario of the RH neutrino dark matter in the context of the minimal 
gauged U{1)b-l model. We have introduced a discrete Z 2 parity in the model, so that one 
RH neutrino assigned as Z2-odd can be stable and, hence, the DM candidate, while the other 
two RH neutrinos account for neutrino masses and mixings through the seesaw mechanism. 
No additional degrees of freedom are necessary to be added. We have evaluated the relic 
density of the dark matter particle. The dominant annihilation modes are via the Higgs 
boson exchange processes in the s-channel and thus, our model can be called Higgs portal 
DM model. It has been found that the relic density consistent with the current observation 
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can be achieved only when the annihilation processes are enhanced by Higgs resonances. 
Therefore, the mass of the RH neutrino DM should be around a half of Higgs boson masses. 
We have also calculated the elastic scattering cross section between the DM particle and a 
proton and found it within the reach of future experiments for the direct DM search. 

Appendix A: The Higgs sector 

The Higgs potential @ contains five parameters: m\ , m 2 , Ai, A2 and A3. These parameters 
can be rewritten in terms of two Higgs VEVs, two physical Higgs masses and the mixing 
angle between them. The stationary conditions are 



ml + Xiv 2 + ^A 3 i/ 2 = 0, 
m\ + \ 2 v 2 + ^A 3 */ 2 = 0. 



(Al) 
(A2) 



The physical Higgs masses are given by Eqs. (jHJ) and (EJ) with the mixing angle that 9 satisfies 

A3W 



tan2 * = -(W-A 2 ^y 

Higgs self interaction terms are expressed as 



(A3) 



C int = \iv^ + \ 2 v'^ + h 3 (v^ 2 + v'^ 2 ) + i(A 1( /) 4 + A 2 ^ 4 + A 3 0V), (A4) 
in terms of </> and i/j. With Eq. (jTJ), these are rewritten in terms of h and H with 9 as 



Hi 



int 



\\v cos 3 9 — A 2 t> ' sin 3 9 + -X^{v cos 9 sin 2 9 — v sin 9 cos 2 9) 



hhh 



+ 



3Aiw cos 2 9 sin 9 + 3A 2 t/ sin 2 9 cos 9 + - A 3 (>(sin 3 9-2 cos 2 9 sin 9) 



+t/(cos 3 9-2 sin 2 9 cos 9))] hhH 
+ 



1 



?>\\v cos 9 sin 2 9 — 3A 2 t>' sin 9 cos 2 + - A 3 (t>(cos 3 9 — 2 sin 2 6* cos 9) 



+v'(- sin 3 + 2 sin cos 2 0))] hHH 



+ 



Aif sin 3 6* + A 2 f ' cos 3 9 H — A 3 (t> sin 9 cos 2 6 + v' sin 2 # cos 1 

2 



+four point interactions. 



(A5) 



We can read off a Higgs three point vertex from Eq. flA5j) . 



In the expression of annihilation cross section, we used the following notations 



d<S> 1 

"77T~ — — 7= COS 6, 

dh ^2 
9$ 1 

m 1 . „ 

^ 1 , A „v 

m=v = 2 cos9 - (A6) 

Appendix B: Amplitude 

We give explicit formulas of the invariant amplitude squared for the pair annihilation 
processes of the RH neutrinos. 



1. Annihilation into charged fermions 



\M 
32 



9b-lWn 



s - Mf, +iM z/ r z > 



+16A 2 



N 



Vf 



<9$ 



(s - Am N 



2 V2 1 



+ 



<9$ 



1 

2 V4 



cos 2 6 



dhs-M 2 + iM h T h dh dH s-M 2 H + iM H T H dH 



(s - Am 2 N ) - mf) . 



(Bl) 



2. Annihilation into neutrinos 



a. Annihilation into v a ^v a (light active-like neutrinos) 



\M 
32 



9b-lWn 



s - M 2 + iM Z 'Y zl 



Am 2 N ) 



3 

8 S 



2U + m - 



)+!(£ + <) cos 2 ^.(B2) 
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b. Annihilation into v s ,v s (heavy sterile-like neutrinos) 



32 



9b-lWn 



s - Ml, +iM z >T z , 



Z' 

dm 



0# <9^ 



dhs-M 2 + iM h T h dh dH s - M 2 H + iM^r H dH 



(s-Am 2 N )(s-Aml s ). 

(B3) 



3. Annihilation into W + W 



\Aif 



8A^ ( ^g 2 v 



dm 



d<p ^ <M> 



d<j) 



dhs-Ml + iM h Yh dh + dH s- M 2 H + iM H Y H dH 
1 



< S - 4 ^»( 1+ 2^(|-<') 2 ) 



(B4) 



Annihilation into ZZ 



\M\ 2 = 8X 2 N (\(g 2 + g' 2 )v 

(s-H) (i + 



90 <9# 
+ 



a0 



i 



2M| V2 



dh s - M 2 + iM h T h dh dH s - M 2 H + iM H T H dH 

M 2 Y) . (B5) 



Annihilation into hh 



Aii denotes the amplitude by s-channel Higgs bosons h and H exchange, while Ai2 does 
that for t(w)-channel N exchange diagram. The formulas for NN — > hH and HH can be 
obtained by appropriate replacement of the vertexes, e.g., Xhhh —> XhhH- 



\Ai\ 2 
\Aii\ 2 



IM1 + M2I 2 , 

dm i 



dh s-M 2 + iM h T h 



iX 



hhh 



+ 



dm 



dH s 



M 2 H + iM H T H 



iX 



Hhh 
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dcos9 2 4 / 
\M 2 \ - A^l — 



/ 



/ 2 

dcos 9 



MiM* 2 = Am N X 3 N 



-8 - I 2 2 + J22 In 



A + 2b 



A -2b 



dh 



dhs-M 2 + iM h T h 



i^hhh + 



-4 + 



s - Am 2 N + A 
2b 



In 



h 

A + 2b 



dHs-Mf I + iM H T H 



(B8) 
i^Hhh 



A -2b 



(B9) 



where 9 is the scattering angle in the center of mass frame. The auxiliary functions appear 
above are defined as 



J22(s,m h ) 

A(s,m h ) 
b(s,m N ,m h ) 



(A + 2a) 2 - 2(s + Am 2 N )A - s(A + m 2 N ) - 3m 2 N (s - Am 



N) 



A 2 - Ab 2 



(BIO) 



— (2A(A + 2a) - A(s + Am 2 N ) + A 2 - Aa 2 - (s - 2m 2 N )(m 2 N - m\) 
+3m 2 N (s - Am 2 N )) , (Bll) 



2 + m h, 



s 2 / ^ 2 



(B12) 
(B13) 



Appendix C: Thermal averaged annihilation cross section 



In partial wave expansion, the thermal averaged cross section is given by 

1 



(av) 



rn 



N 

6 — 

ds 



3 / , . o du> \ T 
w(s) — - I 2w(s) — Am N —r- 



2 
T 



ds J m N 



s=4m 2 N 



s=4m% m N 



with 



Aw{s) = J dLIPS \ M \ 2 = ^ 
where mfi na i is the mass of final state particle. 
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